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摘要 : HSB (Taca) 植物 种 间 在 形态 上 差别 不 大 ， 导 致 分 类 上 存在 一 定 的 困难 。DNA 条 形 码 是 一 种 
利用 短 的 DNA 标准 片段 来 鉴别 和 发 现 物种 的 方法 。 本 研究 利用 核 基因 ITS 片段 和 叶绿体 基因 irnH-psbA， 
rbcL, matK 片段 对 药 昔 暮 属 6 个 种 的 DNA 条 形 码 进行 研究 ， 对 4 个 DNA 片段 可 用 性 ,种 内 种 间 变 异 ， 
barcode gap 进行 了 分 析 ， 采 用 Tree-based 和 BBA 两 种 方法 比较 不 同 序 列 的 鉴定 能 力 。 结 果 显 示 : 单 片 段 
ITS 正确 鉴定 率 最 高 ， 片 段 组 合 rbcL+matK IEW ER, xd CBOL 植物 工作 组 推荐 的 条 码 组 合 rbcL+ 
matK 可 作为 葛 昔 昔 属 物种 鉴定 的 标准 条 码 ， 建 议 ITS FREMRI., mA Tacca integrifolia X 
西藏 的 居 群 与 马来西亚 居 群 形成 了 2 个 不 同 的 遗传 分 文 ， 且 两 者 在 形态 上 也 存在 一 定 的 差异 ， 很 可 能 是 
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一 个 新 种 。 
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Using DNA Barcoding in Genus Tacca (Dioscoreaceae) 
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Abstract: Tacca ( Dioscoreaceae) is a small genus which is taxonomically difficult because of a lack of obvious 
morphologically differences among species. DNA barcoding is a new method of rapid species identification and dis- 
covery using short, standardized genes or DNA regions. In this study, four candidate DNA barcoding regions that 
three (rbcL, matK, trnH-psbA) from chloroplast genome and one (ITS) from the nuclear genome, were evaluated 
among 36 accessions representing 6 species of Tacca. The capability of those regions were evaluated by intra- and in- 
ter-specific divergence and barcode gap analyses, and the identification efficiency was assessed using Tree-based and 
BBA methods. The results indicated that both ITS and the core barcode rbcL+matK proposed by CBOL exhibited the 
highest resolution as a single and combined data respectively. Based on overall performance; matK+rbcL can be 
considered as a potential barcode for identifying the species of Tacca, ITS can be used as a supplementary barcode. 
DNA barcoding revealed two distinct lineages of T. integrifolia distributed allopatrically in Tibet and Malaysia. And 
these two lineages with morphological variations may potentially represent new species. 
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vy 35 SR (Tacca) JRF SERE ( Dioscoreac- 分 布 型 ， 分 布 于 印度 - 马 来 地 区 及 中 南 半 岛 ， 向 
eae) ， 全 世界 约 有 10 种， 中 国有 6 种 (Drenth, 南 延 至 印度 尼 西 亚 诸 品 ， 向 北 延 至 我 国 西 南 、 华 
1972; 335, 2006), BIB BARA, Jez tip 南 地 区 , 个 别 种 分 布 到 澳洲 、 非 洲 (T. leontope- 
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taloides) ， 或 特产 于 南美 热带 (T. parkeri) (Drenth, 
1972) 。 药 碳 昔 属 植 物 花 的 形态 奇特 ， 大 和 苞 片 包 
围 着 聚 伞 状 的 伞 形 花序 ， 并 有 许多 长 须 状 的 小 苞 
片 ， 这 些 特 征 使 得 它们 在 花 开 园艺 贸易 中 广 受 青 
睐 。 药 昔 暮 属 植物 根茎 和 信 药 ， 性 凉 味 若 ， 能 清热 
解毒 ， 消 炎 止 痛 ， 主 治 胃 及 十 二 指 肠 溃 疡 、 高 血 
k, FK., AW, AHR, BAT, RAWE 
的 药 用 价值 ( 贾 敏 如 和 李 星 炜 , 2005 ) 。 由 于 药 
昔 蓝 属 植物 潜在 的 药 用 价值 及 生境 破碎 化 、 片 段 
化 ， 使 得 该 类 植物 居 群 数量 日 渐 缩 小 ， 许 多 种 已 
近 濒 危 ， 大 都 被 列 为 稀有 保护 植物 ( 傅 立 国 和 
金 鉴 明 , 1992) , 

药 昔 车 属 植物 的 系统 学 位 置 和 分 类 标准 一 直 
争议 很 大 (Drenth, 1972; Grierson 和 Green, 1996; 
凌 萍 萍 , 1985), ， 直 到 2003 年 APG II 分 类 系统 将 
本 属 置 于 蓝 蔬 科 中 才 确 定 了 其 系统 位 置 (APG 
II, 2003) 。 由 于 该 属 植 物 分 布 范围 广 ， 不 同 地 
区 形态 上 有 很 大 差异 ， 经 常 造成 该 属 分 类 学 上 的 
混乱 。 另 外 ， 如 老虎 须 (T. chantrieri) , JAEL 
354 (T. subflabellata) 3c 38383 (T. ampli- 
placenta) 等 在 营养 体 阶段 也 很 难 区 分 ， 因 此 果 
实 、 气 孔 、 染 色 体 等 曾 一 度 作为 标准 进行 分 类 
( 姿 萍 萍 , 1985) 。 鉴 于 形态 鉴定 在 分 类 上 的 局 限 
HE (Knowlton, 1993; Jarman 和 Elliott, 2000) , 
目前 迫切 需要 一 种 既 快 速 又 能 弥补 形态 鉴定 的 缺 
陷 的 方法 。 

2002 年 ，Tautz 等 (2002) 首先 提出 将 DNA 
序列 用 于 分 类 学 研究 ， 即 利用 一 段 标准 的 DNA 
序列 作为 标记 来 实现 对 物种 的 快速 准确 的 鉴定 ， 
随后 ，Hebert 等 (2003 ) 正式 提出 了 DNA 条 形 
码 的 概念 ， 并 利用 线粒体 的 COL 基因 序列 对 动 
物 界 的 条 形 码 分 类 进行 了 研究 。 由 于 植物 线粒体 
进化 速度 慢 和 核 基因 的 多 拷贝 等 原因 ， 植 物 的 
DNA 条 形 码 多 从 叶绿体 基因 组 中 选择 ( Chase 
等 , 2005; Cowan 等 , 2006) 。 植 物 条 形 码 片段 的 
筛选 受到 各 种 因素 的 制约 ， 例 如 引物 通用 性 ， 扩 
增 和 测序 的 难 易 以 及 某 片段 在 特定 植物 中 可 能 存 
在 缺失 等 ,不同 的 学 者 提出 了 不 同 的 片段 组 合 
(Chase 等 , 2007; Kress 和 Erickson, 2007; Lahaye 
等 , 2008), lll ITS+tmH-psbA (Kress ^5, 2005) , 
matK + atpF-atpH + psbK-psbl 和 matK + atpF-atpH + 
trnH-psbA 两 个 组 合 (Kim unpublished, 2007), 















































CBOL 植物 工作 组 (2009) 提出 了 rbcL+matK F 
段 组 合作 为 陆地 植物 的 核心 条 形 码 。 

AKA WASH 35 ES B 6 个 种 的 36 个 个 体 为 
材料 ， 对 四 个 DNA Z&JÉ RS Fr Bz ITS, matK, rbcL 
和 trn H-psb A 的 效用 进行 评价 ， 并 尝试 筛选 出 适 
用 于 该 属 植物 物种 鉴定 的 的 DNA. 条 形 码 片段 或 
片段 组 合 ， 为 该 属 植物 的 数字 化 鉴定 及 物种 保护 
提供 参考 依据 。 





1 材料 与 方法 
1.1 研究 材料 

本 研究 选取 了 药 昔 蔓 属 6 个 种 的 36 个 个 体 为 材料 ， 
每 个 种 包括 了 2 ~8 个 个 体 ， 材 料 信息 见 表 1。 实 验 材料 
均 为 野外 采集 后 经 硅胶 快速 干燥 的 叶片 以 备 DNA 的 提 
取 。 赁 证 标本 分 别 存放 于 中 国 科 学 院 西 双 版 纳 热带 植物 
园 标本 馆 (HITBC), ， 中 国 科学 院 昆 明 植物 研究 所 标本 
馆 (KUN) 和 美国 国家 标本 馆 (NMNH ) 。 

12 ”研究 方法 

DNA 提取 采用 CTAB 法 提取 全 基因 组 DNA ( Doyle 
All Doyle, 1987) , 将 DNA F TE 溶液 中 ， 终 浓度 大 约 
J 50 ng: uL', PCR 扩 增 反应 在 PERKIN ELMER (PE) 
9700 型 PCR 仪 上 进行 。 本 研究 中 各 DNA 片段 使 用 的 
PCR 扩 增 和 测序 引物 相同 (3€ 2), matK, trnH-psbA 和 
rbcL 三 个 叶绿体 片段 扩 增 反应 程序 均 采用 : 97%C 预 变 性 
3 min, 94C 变性 1 min，5$2% 退火 1 min, 72°C 延伸 1 
min, 32 个 循环 后 ，72% 延伸 10 min; ITS 片段 扩 增 反应 
程序 采用 : 97°C 预 变性 3 min, 95°C 变性 2 min, 54% 退 
KAS min, 72°C 延伸 1 min, 32 个 循环 后 ，72% 延伸 7 
min。 扩 增 反应 采用 20 pL 的 反应 体系 : 包括 30 ~ 60 ng 
的 模板 DNA, 2 pL 10xbuffer，0.5 单位 的 Taq DNA 聚合 
酶 (5 U-ph'), 1.5 uL ÉJ MgCl, (25 mmol:I)，1.5 
uL ff] dNTPs (10 mmol: L* )， 正 反 向 引物 各 2 pL (5 
pmol:L'!), ， 再 添加 蒸馏 水 至 20 pL。 扩 增 得 到 的 PCR 
产物 经 琼脂 糖 凝 胶 电 泳 检测 合格 后 ， 用 上 海 生物 工程 技 
术 服 务 有 限 公 司 的 纯化 试剂 盒 纯 化 并 测序 。 为 保证 序列 
的 准确 性 ， 所 有 片段 均 采用 双向 测序 。 

HT AUI 33 2B ITS 片段 使 用 引物 ITSA/TTSS 时 
PCR 扩 增 失败 ， 所 以 我 们 也 使 用 了 ITS1 替代 ITS5 用 于 
PCR 扩 增 ( 表 2)。 

1.3 数据 分 析 

用 Contigexpress (version 10.1, Invitrogen) 进行 序列 
的 拼接 ， 拼 接 好 的 序列 用 Clustal X2.0 (Larkin 等 , 2007) 
进行 序列 的 比 对 并 进行 人 工 校正 。 将 4 个 片段 的 DNA 数 
据 进行 单独 及 联合 构建 和 矩阵， 联合 矩阵 中 去 除 序列 不 足 
的 个 体 。 












































































































































表 1 采样 信息 表 


Table 1 Samples for testing potential barcodes 

































































































































































Taxa Locality Voucher Code Taxa Locality Voucher Code 
JK FA-E Tacca plantaginea KRH y T. ampliplacenta 

1 ”广西 桂林 Guilin Guangxi XTBGZL001 T. plantaginea 001 1 PH ALT. Yingiang Yunnan XTBGZL039 T. ampliplacenta 039 

2 广西 桂林 Guilin Guangxi XTBGZL002 T. plantaginea 002 2 RWL Yingjiang Yunnan XTBGZL040 T. ampliplacenta 040 

3 广西 全 州 Quanzhou Guangxi XTBGZL005 T. plantaginea 005 3 ”云南 芒 市 Mangshi Yunnan XTBGZL043 T. ampliplacenta 043 

广西 全 州 Quanzhou Guangxi XTBGZL006 T. plantaginea 006 4 ”云南 潞 西 Luxi Yunnan XTBGZL044 T. ampliplacenta 044 

5 ”贵州 荔 波 Libo Guizhou XTBGZL009 T. plantaginea 009 老虎 须 T. chantrieri 

6 — SUM XE Libo Guizhou XTBGZLO10 T. plantaginea 010 1 AWL Yingjiang Yunnan XTBGZL037 T. chantrieri 037 
ZINIS T. leontopetaloides 2 m AIL Yingjiang Yunnan XTBGZL040 T. chantrieri 040 

1 ”云南 动 腊 Mengla Yunnan XTBGZL003 T. leontopetaloides 003 || 3 广西 宣 州 Xuanzhou Guangxi XTBGZL048 T. chantrieri 048 

2 ”云南 动 腊 Mengla Yunnan XTBGZL004 T. leontopetaloides 004 || 4 云南 普洱 Pu'er Yunnan XTBGZI 049 T. chantrieri 049 

3 ”泰国 清 迈 Chiang mai Thailand XTBGZLO15 T. leontopetaloides 015 || 5 云南 动 腊 Mengla Yunnan XTBGZL050 T. chantrieri 050 

4 ， 泰国 清 迈 Chiang mai Thailand XTBGZLO16 T. leontopetaloides 016 || 6 云南 动 腊 Mengla Yunnan XTBGZL051 T. chantrieri 051 

5 ”马达 加 斯 加 Madagascar XTBGZL017 T. leontopetaloides 017 7 ”云南 马 关 Maguan Yunnan XTBGZL052 T. chantrieri 052 

6 ”马达 加 斯 加 Madagascar XTBGZLO018 T. leontopetaloides 018 8 ”云南 马 关 Maguan Yunnan XTBGZL053 T. chantrieri 053 

7 ”美国 夏威夷 Hawaii USA XTBGZL020 T. leontopetaloides 020 || 2£ 2538 $3 3. T. integrifolia 

8 ”美国 夏威夷 Hawaii USA XTBGZL021 T. leontopetaloides 021 1 马来西亚 芙蓉 市 Seremban Malaysia XTBGZL022 T. integrifolia 022 
Jj erst 38 T. subflabellata 2 马来西亚 芙蓉 市 Seremban Malaysia XTBGZL023 T. integrifolia 023 

1 ”云南 河口 Hekou Yunnan XTBGZLO11 T. subflabellata 011 3 马来西亚 芙蓉 市 Seremban Malaysia XTBGZL024 T. integrifolia 024 

2 ”云南 河口 Hekou Yunnan XTBGZLO12 T. subflabellata 012 4 ”西藏 墨 脱 背 崩 Beibeng Motuo Tibet XTBGZL029 T. integrifolia 029 

3 ”越南 老 街 Laocai Vietnam XTBGZL013 T. subflabellata 013 5 ”西藏 墨 脱 背 朋 Beibeng Motuo Tibet XTBGZL030 T. integrifolia 030 

4 ”越南 老 街 Laocai Vietnam XTBGZL014 T. subflabellata 014 6 ”西藏 墨 脱 汗 密 Hanmi Motuo Tibet XTBGZL007 T. integrifolia 007 
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表 2 四 个 片段 的 引物 信息 表 
Table 2 Primer pairs used for the four DNA markers ITS, 
rbcL, matK, and trnH-psbA 


Region Primers (sequence 5'-3') 





ITS5: GGAAGTAAAAGTCGTAACAAGG 
ITS ITS1; TCCGTAGGTGAACCTGCGG 
ITS4; TCCTCCGCTTATTGATATGC 


1F; ATGTCACCACAAACAGAAAC 


m 724R; TCGCATGTACCTGCAGTAGC 
"T psbA; GTTATGCATGA ACGTAATGCTC 
TRS tmH; CGCGCATGGTGGATTCACAATCC 

390F: CGATCTATTCATTCAATATTTC 
matK 


1326R; TCTAGCACACGAAAGTCGAAGT 


衡量 DNA 条 形 码 的 一 个 重要 指标 是 “barcode gap" 
是 否 存在 ， 即 序列 的 种 内 变异 要 足够 小 而 种 间 变 异 要 足 
够 大 ， 从 而 在 中 间 形 成 一 个 明显 的 间隔 区 (Meyer 和 
Paulay, 2005) 。 本 文 用 Meier 等 (2006) 的 TaxonDNA 软 
件 结合 一 般 的 统计 软件 构建 出 能 够 体现 种 间 、 种 内 的 遗 
传 距离 的 分 布 频 度 的 柱 形 图 。 通 过 Wilcoxon 秩 和 检验 来 
检验 序列 之 间 的 变异 性 差异 。 

用 MEGA 4.0 (Tamura 等 , 2007) 软件 计算 各 DNA 
片段 的 种 内 和 种 间 的 Kimura-2-parameter ( K2P) 成 对 遗 
传 距离 。 同 时 构建 单独 片段 以 及 片段 组 合 的 邻接 树 
(NJ), WE {ELLA 5 000 次 重复 运算 。 参 考 Hollingsworth 
等 (2009) conspecific individuals ‘ grouping together’ 原 
则 ， 即 每 个 种 的 所 有 个 体 都 聚 为 一 个 单 系 分 支 则 被 认为 
该 种 鉴定 正确 。 另 外 ， 利 用 TaxonDNA (Meier 等 , 2006) 
软件 中 的 “Best match”、“Best close match” 和 “All spe- 
cies barcodes” (简称 BBA 法 ) 对 物种 的 鉴定 率 进 行 评估 。 

KRRP, HFW (T. integrifolia) 的 样 
品 在 所 有 的 DNA 片段 中 均 聚 为 2 AS IBIBUA ARA X, H 
了 避免 与 传统 分 类 混淆 ， 我 们 在 本 文中 将 两 个 分 支 处 理 
为 丝 须 药 昔 暮 马 来 西亚 分 支 和 丝 须 药 昔 昔 西 藏 分 文 ， 连 
同 其 余 5 个 种 BIE EA, KRW, KEHE, Ap 
MSS, PALM SSS) 对 条 码 进行 评价 。 







































































2 结果 
2.1 PCR 扩 增 和 测序 成 功率 

PCR 扩 增 和 测序 成 功率 存在 较 大 的 差异 ， 
trnH-psbA 的 扩 增 效率 最 高 ， 为 100% ，rbcL 和 
matK 的 PCR 扩 增 成 功率 相同 ， 均 为 94% ，ITS 
的 PCR 扩 增 成 功率 最 低 ， 为 78% ( 表 3)。PCR 
扩 增 成 功 的 样品 其 测序 率 均 为 100% 。irnH- 
psbA, rbcL 和 matK 三 个 叶绿体 片段 所 有 种 中 均 
能 扩 增 成 功 ， 只 是 个 别 个 体 扩 增 失败 。ITS 片段 
PCR 扩 增 成 功率 低 的 主要 原因 是 裂 叶 欧 昔 暮 所 
有 样品 均 扩 增 失 败 ， 尽 管 我 们 使 用 了 两 对 ITS 引 
物 (ITSSZITS4 fil ITS1/ITS4) ， 并 对 PCR 反应 条 
件 进行 了 优化 。 
2.2 序列 信息 

所 选用 的 4 个 DNA 片段 的 长 度 、 变 异 位 点 
等 信息 见 表 3, ITS 片段 的 矩阵 长 度 包 括 740 个 
位 点 ， 其 中 含有 181 个 变异 位 点 (24.5%) 和 6 
个 插入 /缺失 位 点 ， 老 虎 须 所 有 样品 在 106 ~ 110 
位 点 处 的 有 一 个 4bp 的 插入 /缺失 ; irnH-psbA Fr 
段 的 矩阵 则 包括 349 个 位 点 ， 其 中 22 个 为 变异 
位 点 ， 还 包括 4 个 插入 /缺失 位 点 ， 而 这 4 个 揪 
入 /缺失 位 点 均 为 种 特有 位 点 ， 其 中 丝 须 东欧 警 
的 马来西亚 居 群 在 54 ~ 61 位 点 含有 一 个 8 bp 的 
缺失 ， 裂 叶 药 昔 昔 在 117 ~ 126 位 点 有 一 个 10 bp 
的 插入 ,水 田 七 在 132 ~ 142 位 点 有 一 个 11 bp 
Mika, Bitsy gy EE 144 ~ 150 位 点 含有 一 个 
7 bp 的 插入 ; matK 基因 的 矩阵 长 度 为 742 个 位 
点 ， 包 括 了 62 个 变异 位 点 (8.36%) 和 1 个 6 
bp 插入 /缺失 位 点 ， 这 个 插入 /缺失 位 点 为 丝 须 
东欧 荀 马 来 西亚 居 群 特有 ; rbcL 基因 的 序列 长 
度 在 种 间 无 差异 ， 均 为 651 bp， 其 矩阵 中 仅 含有 
10 个 变异 位 点 (1.54% )。 在 4 个 DNA 片段 中 ， 








表 3 基于 “条 码 分 类 ”的 7 个 种 的 4 个 片段 的 PCR， WF, 变异 位 点 信息 


Table 3 Summary of PCR, sequencing success and variability of the 4 DNA barcode regions. The species resolution, 


intraspecific and interspecific distance was caculated based on the 7 barcode “taxon” names (mean+SE 
F F 





Region PCR Sequencing No. of Indels Aligned Intraspecific Interspecific 
Success/ % Success/% Variable sites/% (length ) length/bp K2P-distance K2P-distance 
trnH-psbA 36(100) 100 22 (6.30) 4 (7, 8, 10, 11) 349 0. 00100. 0004 0.01380. 0028 
rbcL 34(94) 100 10 (1.54) 0 651 0. 0005 x0. 0003 0.0049 x0. 0008 
matK 34(94) 100 62 (8.36) 1 (6) 742 0. 0017+0. 0010 0. 0223 +0. 0027 
ITS 28(78) 100 181 (24.5) 6 (1-4) 740 0. 0033 +0. 0018 0.0903 +0. 0165 
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种 内 遗传 距离 均 明 显 小 于 种 间 遗 传 距离 ;ITS 片 
段 具 有 最 大 的 种 内 遗传 距离 和 种 间 遗 传 距离 ,其 
次 为 matK 和 irnH-psbA ,而 rbcL 则 具有 最 小 的 种 
内 遗传 距离 和 种 间 遗 传 距 离 ( 表 3)。 
2.3 DNA Barcode gap 的 评估 及 种 间 种 内 差异 
的 统计 学 检验 

从 4 个 DNA 片段 种 间 和 种 内 遗传 变异 的 分 布 
图 可 以 看 出 (图 1)，ITS 片段 除了 2.0% ~2.5% 
处 稍 有 重 登 外 ， 种 内 变异 和 种 间 变 异 之 间 的 界限 
明显 ， 表 明 存 在 明显 的 “barcode gap”; 但 是 三 个 
叶绿体 片段 matK , trnH-psbA, rbcL 的 种 内 和 种 内 
遗传 变异 的 分 布 并 不 存在 “barcode gap” (图 1)。 

利用 Wilconxon 秩 和 检验 对 不 同 序列 种 内 种 
间 的 变异 性 进行 统计 分 析 ， 结 果 表 明 ITS 序列 具 
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有 最 大 的 种 间 分 异 率 ， 而 rbcL 具有 最 小 的 种 间 
分 异 率 ， 其 种 间 的 分 异 率 依次 为 ITS > matK > 
trnH-psbA>rbcL ( 表 4); 而 4 个 DNA 片段 的 种 
内 变异 不 存在 显著 性 差异 ( 表 5)。 
2.4 物种 分 辨 率 

本 文采 用 Tree-based 和 BBA 两 种 方法 对 药 莫 
葛 属 的 物种 鉴定 率 进 行 了 分 析 。 基 于 Tree-based 
方法 中 ， 在 基于 单个 片段 的 NJ 树 中 ，ITS 表现 最 
好 (图 2: B), 可 以 区 分 所 有 6 个 种 CERIS] 
昔 暮 扩 增 失 败 没 有 获得 ITS 序列 )，maitK 可 以 区 
分 7 个 种 中 的 5 个 种 (71.43%), trnH-psbA 可 分 
3t7 个 种 中 的 4 个 种 (57.1496) , rbcL 鉴定 能 
最 差 ， 只 能 区 分 7 个 种 中 的 3 个 种 (42. 8696) 















































( 表 6) 。 片 段 组 合 中 除了 rbcL+imH-psbA (57. 1496 ) 
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图 1 4 个 候选 片段 种 间 种 内 K2P 距离 频率 直方 图 . X 轴 表 示 K2P 距离 , Y 轴 表 示 是 出 现 的 累计 频率 


Fig.1 The frequency distribution of inter- and intra-specific Kimura 2-parameter ( K2P) distances for 4 candidate regions. 


The X axis relates to the K2P distances, and the Y axis corresponds to the number of occurrences 












































6 期 赵 月 梅 和 张 玲 : ESS (BRL) 植物 DNA 条 形 码 研究 679 
表 4 4 个 候选 序列 种 间 变异 的 Wilcoxon 秩 和 检验 
Table 4 Wilcoxon signed rank tests of the inter-specific divergence among the 4 candidate regions 
W+ W- Relative Ranks, n, P value Result 
ITS rbcL W=+=20, W-=190, N=20, P<=0. 0007076 ITS>rbcL 
ITS trnH-psbA W=+=36, W-=195, N=21, P<=0. 006035 ITS>trnH-psb A 
ITS matK W+=42, W-=189, N=21, P<=0.01117 ITS»matK 
rbcL matK W+=1, W-=230, N=21, P<=7.424e-05 matK>rbcL 
rbcL trnH-psbA W+=208, W-=23, N=21, P<=0.001385 trnH-psbA>rbcL 
matK trnH-psbA W+=10, W-=221, N=21, P<=0. 0002628 matK»trnH-psb A 
表 5 4 个 候选 序列 种 内 变异 的 Wilcoxon 秩 和 检验 
Table 5 Wilcoxon signed rank tests of the intra-specific divergence among the 4 candidate regions 
W+ W- Relative Ranks, n, P value Result 
ITS rbcL W+=15, W-=0, N=5, P<=0.0625 ITS=rbcL 
ITS trnH-psbA W+=18, W-=3, N=6, P<=0. 1562 ITS = trnH-psbA 
ITS matK W+=12, W-=3, N=5, P«z0.3125 ITS=matk 
rbcL matK W+=1, W-=5, N=3, P<=0.5 matK =rbcL 
rbcL trnH-psbA W+=0, W-=15, N=5, P<=0.0625 trnH-psbA =rbcL 
matK trnH-psb A W+=9, W-=6, N=5, P<=0.8125 matK = trnH-psb A 
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图 2 基于 K2P 距离 的 叶绿体 matK+rbcL 联合 矩阵 (A) 和 核 基 因 ITS (B) 得 到 的 无 根 NJ 树 (分 支 上 显示 高 于 50% 支持 率 ) 
Fig.2  Unrooted neighbour-joining (NJ) tree based on the K2P-distance of cpDNA matK+rbcL (A) 

















and nrDNA ITS (B) used, Bootstrap values over 50% were shown above the branches 
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表 6 基于 “Tree-based” 和 “best match”, “best close match”, “all species barcodes” 的 鉴定 成 功率 
Table 6 Identification successes (% ) based on “Tree-based” and “best match" , “best close match" , “all species barcodes” 
BBA 
Region Tree-based Best match Best close match All species barcodes 
Ps SI AI MI SI AI MI NM SI AI MI NM 

ITS 100 100 0 0 100 0 0 0 100 0 0 0 
rbcL 42.86 41.17 58.82 0 41.17 58.82 0 0 100 0 0 0 
matK 71.43 79.41 20.85 0 79.41 20.85 0 0 88.23 11.76 0 0 
trnH-psb A 57.14 63.88 36.11 0 63.88 36.11 0 0 91.66 8.33 0 0 
ITS+rbcL 100 100 0 0 100 0 0 0 92.50 7.4 0 0 
ITS+matk 100 100 0 0 100 0 0 0 92.59 7.4 0 0 
ITS+érnH-psbA 100 100 0 0 100 0 0 0 100 0 0 0 
matK+rbcL 100 100 0 0 100 0 0 0 93.93 6.06 0 0 
matK-trnH-psb A 71.43 100 0 0 100 0 0 0 94.11 5.88 0 0 
rbcL+trnH-psbA 57.14 100 0 0 100 0 0 0 100 0 0 0 
ITS+rbcL+matk 100 100 0 0 100 0 0 0 92.59 7.4 0 0 
ITS+rbcL+trnH-psbA 100 100 0 0 100 0 0 0 100 0 0 0 
ITS+matK+trnH-psbA 100 100 0 0 100 0 0 0 92.59 7.4 0 0 
matK+irnH-psbA+rbcL 100 100 0 0 100 0 0 0 93.93 6.06 0 0 
matK +trnH-psbA +rbcL+ ITS 100 100 0 0 100 0 0 0 92.59 7.4 0 0 








TE: 成 功 鉴定 (SD; 模糊 鉴定 (AD; 鉴定 失败 (MD ; 没有 匹配 (NM) ; 物种 作为 单 系 的 鉴定 成 功率 (Ps) 





Note; SI; successful identification; AI; ambiguous identification; MI; misidentification; NM: no match; Ps; percentages of species resolved as 


unique clusters 


和 matK+trnH-psbA (71.43% ) 外 ， 其 余 所 有 条 
码 组 合 的 鉴定 率 均 为 100% ( 表 6); BBA 检验 
法 中 ,单独 片段 在 “best match” 和 “best close 
match" 模式 中 表现 为 ITS 最 高 ， 为 100% ， 即 所 
有 序列 均 可 鉴定 到 种 ，matK 为 79. 4196, trnH- 
psbA 为 63. 88% , rbcL 最 小 ， 为 41. 17% ， 所 有 
片段 组 合 均 为 100% 其 中 matK +trnH-psbA 和 
rbcL+trnH-psbA 结果 比 基 于 Tree-based 鉴定 率 更 
高 ， 在 “all species barcodes” 模 型 中 ， 单 独 片 
Pt matK 最 小 ,为 88. 23%， 片段 组 合 中 ITS + 
trnH-psbA, rbcL4trnH-psbA , ITS+rbcL+irnH-psbA 
为 100% ， 其 他 组 合 都 在 90% 以 上 〈 表 6) 。 


3 讨论 
3.1 DNA 条 形 码 的 评价 

DNA 片段 引物 通用 性 是 条 形 人 码 评价 的 重要 
标准 之 一 (Kress 等 ，2005; Hollingsworth ^, 
2009; CBOL Plant Working Group，2009)。 本 研 
究 选 用 的 4 个 DNA 片段 中 ，trnH-psbA 引物 的 通 
用 性 (PCR 扩 增 和 测序 成 功率 ) 表现 最 好 ， 为 
100% ; rbcL 和 matk 各 有 少量 样品 扩 增 失败 ， 
PCR 扩 增 效率 均 为 94% ， 这 可 能 是 因为 样品 的 





DNA 模板 的 质量 较 差 或 实验 技术 方面 的 问题 。 
ITS 的 PCR 扩 增 通用 性 最 低 ， 为 78% ， 主 要 是 
D RUE dE 8 个 样品 中 均 扩 增 失 败 ， 这 可 能 
是 因为 该 种 在 引物 区 核 背 酸 发 生 窗 变 所 导致 
(Liu 等 , 2011) 。 所 有 PCR 扩 增 成 功 的 样品 均 能 
成 功 测序 。 

ITS 序列 具有 最 大 的 种 间 分 异 率 ，maik 次 
Z, rbcL 具有 最 小 的 种 间 分 异 率 ; 而 4 个 DNA 
片段 的 种 内 变异 无 车 性 差异 。 通 过 对 种 内 和 种 间 
遗传 变异 分 布 的 分 析 ，ITS 具有 明显 的 “barcode 
gap”， 物 种 分 辨 率 为 100% ， 而 叶绿体 片段 
matK , trnH-psbA 和 rbcL 的 “barcode gap” 不 明 
显 ， 除 matK Yb, trnH-psbA 和 rbcL 的 物种 分 辨 
率 均 较 低 。 因 此 ， 在 物种 分 辩 率 方面 ITS 片段 作 
为 DNA 条 码 具有 明显 的 优势 ， 但 其 较 低 的 引物 
通用 性 可 能 是 其 作为 该 属 DNA 条 形 码 的 主要 限 
制 因素 ; maik 的 变异 虽然 小 于 ITS, 但 其 高 的 引 
物 通用 性 和 较 高 的 物种 鉴别 率 可 作为 该 属 的 
DNA 条 码 或 条 人 码 组 合 的 一 部 分 。trnH-psbA 序列 
长 度 变异 较 大 ， 含 有 多 个 插入 /缺失 位 点 ， 可 能 
会 对 序列 比 对 造成 困难 ( Hollingsworth 等 ， 
2009) ， 但 利用 好 插入 /缺失 位 点 会 帮助 对 物种 的 
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鉴定 (Kress 和 Erickson, 2007; Liu 等 ，2011 ) ， 
irnH-psbA 在 单独 或 片段 组 合 中 较 低 的 物种 分 辩 
率 ， 可 能 不 适合 作为 该 属 的 DNA RB, KEH 
通用 性 最 高 。rbcL 序列 在 4 个 DNA 片段 中 变异 
最 少 ， 物 种 鉴别 率 也 最 低 ， 但 该 片段 通用 性 较 高 
且 容 易 排序 ， 并 且 在 部 分 片段 组 合 中 具有 较 高 的 
物种 分 辨 率 ， 因 此 可 作为 条 码 组 合 中 的 一 个 候选 
条 码 。 
3.2 药 葛 慕 属 的 物种 分 辩 率 

Tree-based 方法 中 ， 单 片段 鉴定 能 力 为 ITS» 
matK»trnH-psbA»rbcL, ITS 可 以 鉴别 所 有 的 6 个 
种 ， 而 rbcL 的 物种 鉴别 率 为 42. 86% 。 在 片段 组 
合 中 ， 除 rbcL+trnH-psbA 和 matK +trnH-psbA 组 
合 分 别 能 分 辨 57. 14% 和 71. 43% 的 物种 外 ， 其 
它 条 码 组 合 的 物种 分 辨 率 均 为 100% 。CBOL 植 
物 工作 组 (2009 ) 推荐 matK +rbcL 组 合作 为 陆 
地 植物 的 核心 条 码 ， 物 种 分 辨 率 可 达 72% , dE 
本 人 研究 中 ，matK+trbcL 2H & nf AKIA XC CP S 95 
SRA, VIRRAE EXE 10096, MENZ 
种 植物 的 鉴定 的 标准 条 码 。 虽 然 在 所 有 含有 ITS 
片段 的 组 合 条 码 的 物种 分 辩 率 也 均 为 100% ， 但 
由 于 ITS 在 裂 叶 药 昔 薯 中 不 能 扩 增 ， 导 致 该 种 在 
使 用 ITS 时 不 能 有 效 鉴定 。 如 果 能 解决 ITS EA 
Hay SSP PCR 扩 增 的 技术 瓶颈 ,， 那 ITS 单个 
片段 或 与 其 它 条 码 (如 matK All rbcL 等 ) 的 组 合 
也 可 以 作为 该 属 物种 的 标准 条 码 。 由 于 matK+ 
rbcL 具有 三 条 人 码 组 合 和 四 条 码 组 合 相 同 的 物种 鉴 
定 效率 ， 因 此 ， 我 们 推荐 matK+rbel AA E2938 
葛 警 属 植物 鉴定 的 标准 条 码 。 
3.3 ” 丝 须 药草 曹 新 种 的 发 现 

DNA 条 形 码 技术 是 形态 分 类 的 重要 辅助 工 
具 之 一 ， 除 帮助 鉴定 现 有 的 物种 外 ， 还 能 发 现 
一 些 新 的 物种 或 隐 种 (Newmaster 和 Ragupathy , 
2009; Valentini 等 , 2009) 。 由 于 药 莫 暮 属 植物 同 
一 种 的 不 同居 群 在 各 自 特 殊 的 生境 和 气候 地 理 环 
境 的 影响 下 ， 有 些 形态 或 遗传 结构 已 经 发 生 了 变 
化 ,在 物种 的 形态 鉴定 上 缺乏 稳定 的 性 状 〈 张 
玲 等 , 2006) ， 因 此 ， 仅 依靠 形态 特征 来 鉴定 物 
种 具有 一 定 的 困难 ， 需 要 用 DNA 条 形 码 的 方法 
结合 形态 特征 来 重新 评估 现 有 物种 的 分 类 地 位 。 

本 研究 中 ， 基 于 形态 性 状 承 认 的 老虎 须 ， 实 
RAe, RAL gS, TKI, Senay gg 5 























个 种 各 自 均 能 形成 单 系 分 支 ， 得 到 了 DNA 条 形 
码 很 好 的 区 分 (图 2: A, B), mi£Zistg 353 h9 
样品 并 没有 形成 单 系 分 支 ， 而 是 被 分 成 两 个 独立 
的 单 系 分 支 ， 即 采 自 马来西亚 居 群 的 丝 须 药 昔 蓝 
马 来 西 分 支 和 采 自 西藏 居 群 的 丝 须 药 昔 茵 西藏 分 
x (图 2) 。ISSR 分 析 结 果 表 明 ， 丝 须 药 昔 暮 的 
西藏 居 群 和 马来西亚 居 群 间 存 在 显著 的 遗传 分 化 
( 张 玲 等 , 2006) ， 它 们 可 能 已 分 化 形成 了 两 个 不 
同 的 遗传 实体 。 从 分 布 范围 来 看 ， 丝 须 基 昔 莫 分 
布 的 范围 较 广 ， 从 巴基斯坦 、 斯 里 兰 卡 、 印 度 东 
JEN. RSE, Bin, Af, AE, Mee RS 
来 西亚 、 印 度 尼 西亚 群岛 均 有 分 布 ， 且 南北 有 很 
大 的 地 理 隔 离 (Drenth，1972 ) 。 在 我 国 仅 分 布 
在 西藏 东南 部 的 雅鲁藏布江 河谷 内 ， 位 于 热带 雨 
林 的 北部 的 极限 ， 但 由 于 南北 走向 的 河谷 和 其 中 
的 高 山 ， 形 成 了 热带 性 质 的 局 部 小 气候 ( 张 玲 
等 , 2006) ， 这 种 与 中 心 分 布 区 长 期 隔离 势必 会 
造成 遗传 结构 的 分 化 〈Hewitt, 2000) 和 形态 的 
变异 。 丝 须 葛 昔 暮 的 西藏 居 群 和 马来西亚 居 群 在 
形态 上 也 存在 一 定 的 差异 ， 如 西藏 居 群 的 大 苞 片 
颜色 为 白色 ， 下 部 略 紫 ， 而 马来西亚 居 群 的 大 区 
片 为 黑 紫 色 ; 二 者 花期 也 不 同 ， 西 藏 居 群 为 8-9 
月 ， 而 马来西亚 居 群 为 5-6 月 。 综 合 分 析 ， 我 
们 认为 丝 须 药 昔 暮 西藏 居 群 很 可 能 是 一 个 新 种 ， 
但 还 需要 更 详细 的 形态 学 证 据 对 其 分 类 学 地 位 进 
行 评 价 。 本 文 的 结果 对 于 药 昔 暮 属 植物 的 保护 和 
分 类 学 的 认识 具有 重要 的 作用 。 
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